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ABSTRACT 



Most developments of Cerenkov Radiation are in terms of the 
Fourier components of the fields and power emitted by a single 
electron. When many electrons in a compact bunch are emitted from 
an accelerator, the bunch radiates coherently and at a lower 
frequency than for a single electron. The theory for the time 
structure of the fields arising from a charge bunch is developed, 
and it is shown that the source of the radiation is di/dt. 

Present detector technology should be able to resolve these 



fields . 



- 





TIME DEVELOPMENT OF CERENKOV lADI^TION 
INTRODUCTION 

Cerenkov radiation, produced by a charqe or qroup of charqes, 
novinq faster than the speed of electromaqnetic radiation in a 

nedium, has been investiqated , startinq with the exDeriments of 

1 . 2 . 
lerenkov in 1934 and the exDlanation by Frank and Tamm in 1937. 

Since power radiated by a single charoed particle is proportional 

co the frequency, most of the research effort has been devoted to 

the relatively intense optical radiation which is favored over the 

3 4 

nicrowave reqion by a factor of about 10 . The optical results ' 
are qiven in terms of the Fourier components of the fields and the 
radiated power. 

In our previous work^'^ it was noted that microwave radiation 
can be siqnificant because all the electrons in an accelerator 
aunch (about 1 0 9 ) radiate coherently; an effect which more than 
offsets the sinqle particle increase in radiated power with 
frequency. For an electron beam aenerated by a traveling wave 
Linac and passing through air, it was shown that the various 
larmonics of the basic frequency up to about the tenth are emitted. 
[In the case of an L or S band Linac, these correspond to 10 GHz 
ind 30 GHz respectively.) 

The time structure of Cerenkov radiation fields in the optical 
ind even in the higher frequency microwave regions is difficult to 
observe because the detectors reoister power. One of the few 
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treatments of the time dependence, by Tamm in 1939, showed that 
the optical radiation by an electron is singular on the Cerenkov 
front. Here we consider the time structure of fields generated 
when electron bunches radiate coherently; in a development which 
complements the frequency domain analysis of our earlier work 
The fields should be observable for beams from induction 
accelerators which produce bunches much longer than those produced 
by S or L band Linacs. 
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MAGNETIC RADIATION FIELD 



The purpose of this paper is to present a development of the 
time dependence of the electric field generated by the Cerenkov- 
mechanism. The method is first to determine the potentials from 
the moving charge distribution, and subsequently to obtain the 
fields (in cgs units) from the potentials by 



B = V x A 



( 1 ) 



E 




( 2 ) 



We assume a charge density function and a current density 
j v = p^v/c^ with the velocity v m the plus z direction. The 
charge and current are assumed to be concentrated along the z axis 
such that 

P v ( r , t)=p(z,t)5(x)5(y) (3) 



and the charge is assumed to move with no change in shape so that 
the z and t dependence of the charge is 

p ( z , t ) = p (z-vt) (4) 

0 

Note that p^ and 3 represent the usual charge and current 
densities, while p and p Q throughout this paper are charge per 
unit length. The velocity of light is c and c q in the medium and 
free space, respectively. 
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The potentials are found by taking the usual retarded 
solutions to the wave equations; which become under the 
assumption of a line distribution of charge (3), 



4> ( r , 


t) 


= e 1 J R ^(r*, t') dz' 


(5) 


-y . -y 

A ( r , 


t) 


-y 

= / R-ip(r' # t') dz 1 


(6) 



o 



where 



-+■ 

R 



r 



r' and t' 



is the retarded time 



f = t - |r - r ' |/c 

Now (4), the assumption of rigid motion of the charge 
distribution, can be incorporated into the potentials, and a new 
variable u(z') = z' - vt ' can be introduced so that the potential 
(5) and (6) become 



i>(r, t) = e 1 / R (u) dz ' (8) 

J o 

•y 

A (r,t) = -^ / R“ 1 p (u) dzi (9) 

c J o 

0 

/ 

Also, since the charge is confined to the z' axis, the new 
variable u(z') can be written more explicitly. 



u ( z ’ ) = z ’ - 



vt + - [x 2 + 
c 



( z - z ' ) 2 ] 



1/2 



( 10 ) 
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The magnetic field B may be calculated from (1) and since A 
has only a z-component, B has only the x and y components, 

3 8 

B = — A and B = - — A . Carrying out the differentiation 
x 9y z v 9x z j r 

for the x component gives 



B x = b ! I7 R 1 % (u)az ' 

0 2 



(ID 



+ £- 1 R "' I7 % (u) dz ' 

0 2 



For radiation, the first integral, falling off as R“ 2 at large 
distances, will be neglected and only the second term will be 
considered further. From (10), it is seen that u is a function 
of x and y so that the second integral can be written 



B 



x 




p ' ( u ) dz ' 
o 



(12 ) 



where p '(u) is the derivative of p with respect to its 
c c 

argument u. The corresponding expression for B has y replaced 

by (-x). These two components can be combined to give the 

. 

total magnetic radiation field B. In the cylindrical 

coordinates, (s, 9, z) where s is the radius vector 

1 / 2 + 

s = (x 2 + y^) , B is tangential (i.e. in the 6 direction) 

with a magnitude given by 



v z s 

B = / - 

cc R 



p ' ( u ) dz ' 
2 0 



(13 ) 
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TIME DEVELOPMENT 



In order to evaluate (13) for B, it is necessary to consider 

the dependence of u on z' as qiven in (10). in the u-z' plane, 

the first two terms are a straiqht line with unit slooe and an 

intercept which chanaes with time, while the third term is a 

hyperoola openinq in the +u direction with asvmptotic slopes of 

- — . The sum of these two curves is u(z'). In the Cerenkov 
c 

case with v > c, the result is a curve whose ends both point 

upward as shown in Fiq. 1. As time increases, the entire curve 

will translate downward to smaller u values as a result of the 

neqative second term in (10). 

Only chancinq currents (those with a non zero p ' will 

o 

contribute to the maqnetic radiation field (13). To proceed 

and demonstrate the method, a ramo-front current pulse is 

chosen as a simple example. Assuminq that the front end of a 

current pulse increases linearly up to a constant value, the 

derivative p '(u) will be a constant valued square oulse of 
o 

maqnitude p ' as is also shown in Fiq. 1. The correspond inq 
m 

neqative p q '(u) pulse occurrinq at the tail of a current pulse 

is not shown and its effect is considered separately. 

For large neqative times, the u(z') curve (a) is completely 

above the oulse-like non zero portion of the p q '(u) curve so 

that the contribution in (13) to B from p ' (u) is zero and 

o 

therefore, B is zero. As time increases, the u(z') curve 
moves downwar'd until the B Pulse beqins when u(z') is tanqent 
(curve b) to the upper portion of the p q '(u) pulse. The value 
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of the integral in (13) increases as u(z) continues its 
constant downward motion with increasing time until u(z' ) 
becomes tangent (curve 3) with the lower part of the p ' (u) 
pulse. At this time the non zero part of the integral has the 
largest extent -from zj to z 2 - At later times, the integral 
breaks into two regions of the z' axis and if p q '(u) is 
constant, the value of the integral decreases with increasing 
time because the extent of the integral in the two regions 
continues to decrease as a result of the upward turn of u(z'). 

This calculation may be carried further to determine the 
time structure (shape) of the resulting B pulse. Although the 
expression (13) for B can be integrated directly in the case 
where the slope p q '(u) is constant, it is instructive to carry 
out the calculation by developing u in a power series. 

Denoting z ' as the value of z 1 at which u(z') has zero slope, 
the values of z ' , u(z 1 ) and the second derivative are 



z 



l 

0 



s 





— 1/2 



(14) 



u (z ' ) 
0 





(15) 




= 2 A 



(16) 



so that u can be expressed as a power series about the minimum 

u = u(z • ) + A(z' - z ■ ) 2 (17 ) 

o o 

The limits Zj' and z 2 ' can be written in terms of the 

minimum value as z 1 2 = z ' + Az and zi 1 = z ' - Az ' where Az ' 

L o 1 o 

is the value of z ' -z ' such that the difference 

o 

u(z') - u(z ' ) = a, the width of the current derivative pulse 
o 

p Then from (17), a = A(Az') 2 or 

0 

Az' = (a/A) 1/2 (18) 

Using this value, the maximum magnetic radiation field for the 
rising front of the magnetic field pulse is easily evaluated 
from (13) under the assumption that s and R are slowly 
varying to give 



A ' 

0 

Values of B for the rise up to the peak value given above 
are found by the same process but using appropriately smaller 
values of Az ’ . The result is that the integral (and therefore 
B) increases as t 1/2 after the onset of the pulse. After the 
maximum magnetic field is reached, the integral splits into two 
parts. If the expression on the right side of (19) is called 
1(a), the value of B at later times becomes 
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B = I (a ’ + a) - I (a ' ) 



( 20 ) 



where a' is the distance by which the minimum in u(z') is below 
the lower step of the p Q ' pulse in Fig. 1. The first term in 
(20) increases slowly with a', but the second term decreases as 
(a* ) 1/2 leading to the sharp fall off of the magnetic field 
after the maximum as shown in Fig. 2. 

A complete current pulse may be considered as a linear 
rise, followed by a constant current, and then a linear 
decrease. The latter part gives rise to a negative p'(u) and a 
reversed magnetic field pulse so that the magnetic field for a 
complete current pulse has the double peaked structure shown in 
Fig. 2 . 
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ELECTRIC RADIATION FIELD 



In a manner similar to the derivation of (13), the electric 
radiation field may be found from (2), (8), and (9). The 

details are omitted, but the result is 



2 



V 

c 

0 




R c 
R c 




p ’ ( u ) dz ’ 
o 



(21 ) 



The direction of E may be determined from the following 
considerations. If & is assumed approximately constant and 
denoted by R in the region which contributes most strongly to 
the integral, then 



E 





I (1 - —cos 9) 
c 



( 22 ) 



where I represents the integral in (21) without the factor in 
• . 

parenthesis and 9 is the angle between R and the z axis. 

m 

But the value of R in the region which contributes to the 

m ’ 

integral is found by evaluating the general expression (21) at 

z’ = z q ' . To simplify the expression, let the’ observer be at 

z = o and also assume that the p ' pulse is centered near 

u = o. Then R = (s 2 + z ' 2 ) 1 / 2 may be evaluated using (14) to 
m 

gi ve 




( 23 ) 
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If the usual Cerenkov angle is defined as cos 9 = — 

c v 

can be written as 



R 

m 



R 

m 



s 

sin 9 

c 



(24) 



x *> -V 

From (22 J it is apparent that E is perpendicular to R^ when 

9=9 and (24) shows that R is the value of R at the Cerenkov 
c m 

angle 9 . Thus, the electric field from the front of the pulse 
c 

(i.e., z' = z ) is transverse to R . 

m m 

The situation is clarified in Fig. 3. The charge, 

traveling from left to right, emits a signal from A, which 

travels to the the observer at 0, traversing a distance R . 

m 

The observer is at z = o and a distance s from the path. The 

field E is perpendicular to R . The signal was emitted from A 

at an earlier time t' in order to arrive at 0 at the time t, 

with c(t-t') = R . By the time the signal reaches 0, the 

particle is at B, with D = v(t-t' ). Then R /D = c/ v = cos9 as 

m e 

expected. D is the path length from A to B. 

From Fig. 3, one should also note that the field is 

transverse to R , which points from the earlier (retarded time) 

m i 

position of the particle, and is radial relative to the present 

* 

position of the particle. The former condition holds for 
typical dipole radiation, while the latter condition holds for 
the Lienard-Wiechart field for a particle moving with v < c. 

The Lienard-Wiechart fields fall off as the inverse square of 
the distance, and do not represent radiation. In contrast, the 
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fields discussed here fall off more slowly than R “ 1 and 
represent radiation; and the total radiated power is discussed 
in the next section. 
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RADIATED POWER 



The energy radiated may be found bv calculating the 
Povntinq vector and integrating over a surface. If the surface 
is a cylinder centered on the z axis, the fields at a given 
time have a pattern independent of angle and a z dependence as 



shown in 


the top 


curve of 


Fiq . 2 . 


The Povntinq vector is. 


of 


course , 


alonq R , 
m 


and the 


outward 


component may be 


inteqra 


ted 


over the 


cylinder 


to give 


the total power radiated. 


As a 


crude 


estimate 


for the 


integral , 


replace 


the field bv the 


peak f 


ield 


(19) and 


let the 


spatial width of 


the Pulse be a. 


The rad 


iated 


power is 


then 


P ( aporox ) 


R 

c 2 


vl 2 sin 2 9 
o c 




( 25) 



in cqs units. In the mks system, the sauare bracket is 
replaced by 2uo/u. 

In the earlier calculations^'^ the fields and power were 
exDressed in terms of Fourier amplitudes. If the same current 
pulse is assumed, P has frequencv components up to the value 

U) * 

of to such that the wave lenqth of the radiation is eaual to the 

pulse lenqth. If it is assumed that P rises linearlv up to 

to 

this frequency and suddenly drops, the total cower radiated 
becomes (in Mks units) 



( 26 ) 



P =^u vl 2 sin 2 9 
L u) 2 o o c 

Equations (25) and (26) are both rough estimates and the 
point is that the similarity of the results is asserted to be 
confirmation that the calculations in this paper represent the 
Cerenkov radiation, here expressed in terms of time dependence 
of the fields. 
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DISCUSSION 



In precedinq sections, the time structure of the electric 
and magnetic radiation fields was developed. Only the far 
fields were retained and in the development leading to the B 
field (13), and the E field (21), only the assumption of a 
rigid charge confined to a line was introduced. It is seen 
from these equations’ that the time derivative of the current is 
the source function. 

The simple model chosen to demonstrate the method of 

3evelooinq the time structures was that of a uniform charge 

iistribution with uniformly varyinq front and rear sections. 

Chis model qives the square pulse charqe derivative of heiqht 

3 1 shown in Fiq . 1 which is easy to use in evaluating the 

m 

.nteqral (13). Similar remarks hold for the power series 
expansion of u(z') which is an increasingly better 
approximation as the time durinq which the current is chanaino, 
decreases. Current variations other than linear mav be readilv 
incorporated within the framework given. Also it should be 
loted that in all cases the variation of R -2 in (I'M which was 
assumed constant in the example will tend to sharpen both the 
Leading and trailing edges of a field pulse. 

In the evaluation of the time structure of the fields, the 
jeak field arose when the integral (13) had the most widely 
spaced limits; a situation which occurs because u(z') has a 
legative slope for sufficiently negative values of (z' - z) as 



shown at the left side of Fig. 1. In the non-Cerenkov case 

(— < 1) this situation does not arise since then the slope of 
c 

the u(z') function always has the same sign. In this case 
(i.e. v < c) the u(z') curve bends downward instead of upward 
for large negative values of ( z ' - z) and the only contribution 

to the integral (13) is from small regions of z'. 

These results show how the time structure of Cerenkov 
radiation arises from the time rate of change of the charge 
distribution in an electron bunch. Present technology is such 
that this structure is not observable in the Cerenkov radiation 
from S or L band Linacs because of their relatively high 
fundamental frequency. However, induction accelerators with 
their longer electron bunch structure should produce Cerenkov 
signals in air for energies greater than about 25 Mev, which 
should be observable. 

The extension of this method of calculation of the fields 

for both Cerenkov and sub-Cerenkov charge velocities is easily 

made for cases for which the charge derivative p'(u) is not 

constant. A detailed report is under preparation. 

Finally, we note that although the results of other 
4 7 

workers ' often have singularities in the radiated power at 

the Cerenkov angle, the present results and our previous 
5 6 

ones ' show that the radiated power is finite whether 
calculated in the frequency or time domain. Also it should be 
noted that causality is satisfied because the fields are zero 
at times earlier than the leading edge of the pulse shown in 
Fig. 2 . 
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FIGURE CAPTIONS 



Figure 1. The function u = z* - vt ' , defined in the text is 

plotted for increasing times tp 1 2 / at the point 

of observation. The corresponding current derivative 

profile, on the right, is a function of u only and 

remains fixed in time. The field signal pulse starts 

at t , and reaches a maximum at t 3 . 

2 



Figure 2. The electric field pulse generated by the beam 

current profile, shown in the lower curve. 

Figure 3. Geometrical relations for the Cerenkov pulse. The 

source ( p ' ) at A emits a signal at an early time 
0 

giving the 2 field at the observer 9; when the field 
reaches the observer, the particle is at B. 
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Figure 1. 
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